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SUMMARY 


OBJECTIVES 

Examine  the  design  parameters  for  systems  of  multiwire  cages  carrying  high  voltages 
above  ground  in  order  to  minimize  the  maximum  electric  field  on  the  surface  of  the 
wires.  For  the  two-wire  case,  compare  the  following  methods  of  calculation:  (1)  a  simple 
geometrical  theory  formula  (SGT);  (2)  a  complete  geometric  theory  (CGT);  (3)  an  infinite 
series  of  images  (“exact  solution”  for  given  constraints);  and  (4)  a  computer-based 
method  of  moments  (MOM)  solution.  Calculate  the  MOM  results  for  cages  with  2  through 
6  wires,  and  for  12  wires. 

RESULTS 

The  MOM  accurately  matches  exact  solution  results  for  the  two-wire  case.  The  CGT 
formula  reproduces  the  results  of  the  more  exact  solutions  fairly  well.  For  3  or  more 
wires,  a  comparison  of  results  of  MOM,  SGT,  and  CGT  suggests  a  modification  to  the 
SGT  called  modified  simple  geometric  theory  (MSGT).  The  MSGT  requires  few  calcula¬ 
tions,  can  be  applied  to  cages  with  any  number  of  wires,  and  reproduces  the  MOM  results 
surprisingly  well  for  all  the  cases. 

The  results  show  that  there  is  an  optimum  wire  spacing  to  minimize  surface  electric 
field.  This  spacing  depends  upon  the  number  of  wires  in  the  cage,  their  diameter,  and 
height  above  ground.  Tlie  optimum  dimensions  have  been  tabulated  and  design  curves  are 
given  for  cages  using  2,  3,  4,  5,  6,  and  12  wires. 

The  concept  of  electric  field  equivalent  radius  is  found  to  be  extremely  useful.  Appli¬ 
cation  of  this  concept  leads  to  a  simplified  design  procedure.  Two  design  procedures  are 
given,  one  based  on  minimum  wire  spacing  for  a  given  wire  diameter,  and  the  second 
based  on  minimizing  the  total  wire  cross  section  of  the  cage.  These  procedures  are  pre¬ 
sented  with  three  examples  that  include  antenna  and  transmission  line  applications.  Ap¬ 
plication  to  cylindrical  geometry  is  considered  as  is  the  effect  of  shielding  from  other 
objects  in  the  vicinity  of  the  cage. 

The  modified  simple  geometrical  formula  and  the  design  procedures  given  are  useful 
for  design  of  high-voltage  transmission  lines  and  high-power  transmitting  antennas. 
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1.0  INTRODUCTION 


1.1  GENERAL 

The  application  of  multiwire  cages,  sometimes  called  bundles,  to  high-voltage  anten¬ 
nas  and  for  power  distribution  goes  back  many  years.  The  early  workers  in  this  field 
found  that  several  wires  together  could  carry  more  voltage  than  a  single  wire  (Clark, 
1932;  Temoshok,  1948;  Adams,  1955;  Miller,  1956).  In  the  past,  approximate  formulas 
have  been  used  to  calculate  the  surface  electric  fields  on  the  wires  in  a  cage.  By  using 
modem  computer  techniques  to  generate  more  exact  data,  we  have  discovered  how  to 
modify  the  simple  formula  to  provide  more  accurate  results. 

For  a  fixed-wire  diameter,  a  fixed  number  of  wires,  and  a  fixed  cage  height,  there  is 
an  optimum  spacing  having  minimum  surface  gradient. 

1.2  OBJECTIVE 

This  report  will  present  the  modified  formula  for  the  calculation  of  surface  electric 
fields  on  multiwire  cages,  and,  two  optimum  design  techniques  for  multiwire  cages. 

1.3  BACKGROUND 

Many  transmitting  antennas  in  the  lower  frequency  bands  are  of  the  top-loaded 
monopole  type.  Often,  the  top  load  is  composed  of  cables  suspended  from  towers,  or 
across  valleys,  etc.  These  cables  can  carry  high  voltages  and,  in  many  ways,  exhibit  simi¬ 
lar  design  problems  as  the  high-voltage  power  transmission  lines. 

If  the  antenna  is  small  with  respect  to  wavelength,  than  the  maximum  power  it  can 
radiate  is  proportional  to  top-load  area  squared,  frequency  to  the  fourth  power,  and  maxi¬ 
mum  voltage  squared  (Hansen,  1990;  Watt,  1967).  The  maximum  voltage  limit  for  an 
antenna  is  determined  by  voltage  breakdown  on  insulators,  and/or  by  corona  formation  on 
the  structure.  Conductor  selection  for  the  various  portions  of  the  antenna  is  one  of  the 
most  important  design  considerations.  The  conductor  size  for  high-voltage  applications  is 
often  established  by  corona  performance  rather  than  by  current  carrying  capacity  or  struc¬ 
tural  considerations. 

1.3.1  Corona 

Corona  forms  when  the  surface  electric  field  (potential  gradient)  exceeds  the  break¬ 
down  strength  of  air.  Undesirable  effects  of  corona  include  (1)  radio  noise,  (2)  audible 
noise,  and  (3)  power  loss.  In  addition,  chemical  byproducts  of  corona  can  increase  corro¬ 
sion  rates.  Considerable  study  of  these  effects  has  been  done  by  the  power  distribution 
community  (Peek,  1929;  EPRI,  1987).  A  further  consideration  for  high-0  VLF/LF 
antennas  is  that  a  significant  amount  of  corona  will  change  the  effective  antenna 
capacitance  and,  thus,  make  it  difficult  to  maintain  tuning. 


1 


The  voltage  at  which  corona  first  appears  is  called  the  onset  voltage  Vq.  The  larger  the 
wire  size,  the  larger  the  Vo-  Tq  is  considerably  reduced  for  wet  conditions.  Economy 
dictates  that  power  distribution  lines  be  designed  to  operate  with  some  corona,  especially 
during  inclement  weather. 

The  effects  of  corona  are  proportional  to  frequency.  For  VUF,  and  even  more  for  LF, 
having  any  portion  of  the  antenna  in  corona  is  undesirable.  At  these  frequencies,  even  a 
small  segment  of  wire  in  corona  can  dissipate  a  large  amount  of  power  (Watt  &  Hansen, 
1992). 

The  onset  of  corona  correlates  with  the  maximum  surface  electric  field  (Peek,  1929; 
Comber  et  al,,  1987).  For  this  reason,  the  surface  electric  field  (E-field  often  called  volt¬ 
age  gradient  or  gradient)  is  used  to  evaluate  corona  performance.  The  surface  gradient 
resulting  in  corona  onset  is  called  the  critical  onset  gradient  £o-  The  critical  gradient  for  a 
single  conductor  at  60  Hz  has  been  the  subject  of  considerable  study  (Peek,  1929;  EPRI, 
1987).  Peek  developed  an  empirical  formula  for  Eq  versus  wire  diameter  at  60  Hz  for  dry 
single  conductor  wires: 

£o  =  2.107  d{\  +  0.952/(d  a)0-5}  kV/mm  rms, 

where 

d  -  relative  air  density  (1  @  STP),  and 
a  »  wire  radius  (mm). 

This  can  be  applied  to  the  calculation  of  the  onset  of  corona  by  using  a  surface  rough¬ 
ness  factor  M  and  a  frequency  factor  F.  The  surface  roughness  factor  takes  into  account 
that  the  local  gradients  are  higher  than  the  average  gradient. 

Some  study  has  taken  place  of  the  critical  gradient  in  the  VUF  and  LF  bands.  The 
critical  gradient  at  VUF  and  LF  for  dry  wires  follows  the  same  law  as  at  60  Hz,  with  a 
slight  reduction  for  frequency  (Smith,  1963,  1985;  Watt,  1967). 

The  critical  gradient  depends  on  many  variables.  Among  the  important  parameters  are 
frequency,  atmospheric  conditions  (e.g.,  pressure,  temperature,  and  humidity),  surface 
conditions  (wet  or  dry,  oil  or  other  contaminants),  and  surface  roughness.  In  particular, 
for  both  VLF  and  LF,  the  critical  gradient  is  considerably  reduced  for  wet  conditions. 

Wire  diameter  has  a  strong  influence  on  the  critical  gradient,  with  smaller  wires  hav¬ 
ing  larger  critical  gradients,  because  the  electric  field  falls  off  more  rapidly  with  distance 
from  the  surface  for  smaller  wires.  For  ionization  (corona)  to  occur,  electrons  from  the 
surface  must  reach  the  critical  velocity  for  ionization  within  a  few  collision  distances.  This 
critical  velocity  depends  on  the  type  of  molecules  involved  in  die  collisions.  Electron 
velocity  is  proportional  to  the  integral  of  electric  field  times  distance  from  the  surface. 
Because  of  the  rapid  decrease  in  electric  field  near  the  surface  of  smaller  wires,  the 
surface  field  must  be  higher  to  support  corona.  However,  since  the  surface  electric  field  is 
inversely  proportional  to  radius,  the  voltage  required  to  give  this  critical  gradient  is  less 
for  smaller  wires.  Thus,  smaller  wires  go  into  corona  at  lower  voltages,  but  with  larger 
surface  gradients  than  larger  wires. 
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The  variation  of  critical  gradient  with  radius  is  strongly  influenced  by  surface  condi¬ 
tions,  polarity,  and  frequency.  For  example,  large  wires,  when  wet,  have  a  lower  critical 
gradient  than  when  dry,  while  small  wires  can  exhibit  the  opposite  performance. 

The  wire  diameter  used  for  antennas  and  transmission  lines  usually  exceeds  1/2  inch. 
Wires  of  this  diameter  have  a  lower  critical  gradient  when  wet.  The  wet  value  is  used  as 
the  design  parameter  for  outside  applications. 

A  complete  discussion  of  this  phenomenon  is  beyond  the  scope  of  this  paper.  Experi¬ 
ments  have  shown  that  corona  effects  can  be  conveniently  expressed  in  terms  of  the 
maximum  conductor  gradient.  For  this  reason,  the  design  specifications  for  corona  per¬ 
formance-are  stated  in  terms  of  maximum  allowable  surface  gradient.  For  simplicity,  a 
single  design  value  of  surface  gradient  is  often  specified.  This  is  usually  acceptable  for 
the  larger  diameter  wires  typically  used  for  transmission  lines  and  antennas,  because  the 
critical  gradient  versus  diamecer  does  not  change  too  much  over  the  range  of  the  wire 
diameters.  For  VLF/LF  antennas  near  sea  level  with  wires  in  the  vicinity  of  1 -inch-diame¬ 
ter,  common  design  criteria  use  maximum  surface  gradients  in  the  range  of  0.7  to  0.9 
kV/mm  RMS  (Smith,  1985).  For  60  Hz  applications,  a  design  limit  for  wet  conditions 
used  by  some  power  line  designers  is  1.2  kV/mmi  (EPRI,  1987). 

1.3.2  Surface  Gradient  Calculations 

Calculation  of  the  exact  gradient  at  every  point  on  the  surface  of  cables  is  complicated 
by  the  fact  that  the  cables,  often  multiply  stranded,  usually  are  not  smooth.  In  addition, 
the  geometry  generally  causes  the  electric  gradient  to  be  nonuniform  around  the  wire 
periphery. 

To  simplify  the  calculation  of  the  surface  electric  field  for  corona  performance  deter¬ 
mination,  the  conductors  are  usually  treated  as  smooth  cylinders.  Correction  factors  are 
then  used  for  surface  roughness  and  frequency  effects.  For  stranded  wires,  the  typical 
correction  factor  for  surface  roughness  {^f)  varies  between  0.83  and  0.87  (Chaudhari,  see 
Miller.  1956). 

The  effect  of  frequency  on  Eq  is  small  at  VUF  and  LF.  The  frequency  effect  factor  F 
has  been  found  empirically  to  be  on  the  order  of  0.90  at  VLF  (Watt  &  Hansen,  1992). 

To  evaluate  the  corona  performance  of  multiple  wire  cages,  it  is  necessary  to  calculate 
the  value  of  the  surface  electric  field.  In  the  past,  several  authors  have  examined  this 
problem  (Clark,  1932;  Temoshok,  1948;  Adams,  1955;  Miller,  1956;  Comber  et  al.,  1987; 
Doyle  et  al.,  1987).  Miller  (1956)  has  some  measured  and  calculated  60  Hz  results.  Al¬ 
though  some  recent  efforts  to  calculate  the  gradients  on  cages  make  use  of  computer 
techniques,  ail  involve  some  approximations.  For  cases  that  can  be  approximated  by  a 
single  isolated  cage  above  ground,  a  simple  formula  had  previously  been  developed.  In 
the  following  sections,  this  formula,  as  well  as  some  other  formulas,  will  be  developed 
and  compared  with  computer  results  calculated  by  an  “essentially  exact”  method-of- 
moment  technique. 


*  1.  Personal  communication  with  J.  Bainbridge,  Transmission  Line  Engineering  Section,  San  Diego  Gas  and  Electric 
Corporation,  September  1987, 
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1.4  CAGES 


1.4.1  Definitions/Geometry 

The  cage  geometry  is  illustrated  in  figure  1  for  the  genera!  n  wire  cage.  The  cage 
center  height  above  ground  is  (A),  the  radius  of  the  cage  is  (b),  and  the  radius  of  the 
individual  wires  is  (a).  MKS  units  are  used. 


FROM 

GROUND 

Figure  1.  n-wire  cage  geometry. 

1.4.2  Assumptions 

Several  assumptions  are  operative; 

(1)  The  wires  are  parallel  to  each  other  and  to  ground. 

(2)  The  wires  are  perfectly  smooth  circular  cylinders. 

(3)  The  wires  and  the  ground  are  perfect  conductors. 

(4)  Ail  wires  are  at  the  same  potential  V  . 

(5)  The  wires  ail  have  the  same  radius  a. 

(6)  The  wires  are  uniformly  separated  on  a  circle  of  radius  b. 

(7)  The  wires  are  much  farther  above  ground  than  the  cage  radius  b  (i.e.,  h  »  b). 
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The  last  assumption  results  in  some  simplifications: 

(1)  The  charge  is  the  same  on  all  wires. 

(2)  The  electric  fields  in  the  vicinity  of  the  wires  are  not  affected  by  the  presence  of 
the  ground  (note;  charge  magnitude  is  affected). 

1.5  THEORY 


1.5.1  Matrix  Equations 
A.  General  Equations 

The  charge  on  each  wire  is  assumed  to  be  a  filament  located  at  the  center  of  the  wire. 
This  is  a  valid  approximation  if  the  wires  are  separated  far  from  ground  and  from  each 
other  relative  to  their  diameter.  As  the  wires  come  closer  together,  this  approximation 
breaks  down.  This  effect  will  be  examined  in  the  following  section  and  a  correction  term 
developed. 

Following  the  notation  and  method  of  Miller  (1956),  assuming  wires  with  the  charge  at 
the  center,  the  Maxwell  potentials  define  the  relationship  between  voltage  and  charge  as 
follows: 

Self  Potential 

Pii  =  H2h/a;)/{2ji€Q)  , 

Mutual  Potential 

P,j  =  ln(D,y<<„)/(2jte„)  , 

where  a,,  is  the  radius  of  the  /th  wire,  d^j  is  the  distance  between  wires  /  and  j,  while  D^j ,  is 
the  distance  between  wire  i  and  the  image  of  wire  j  reflected  in  the  ground  plane.  Cq  is  the 
permittivity  of  free  space. 

The  relationships  between  voltage  and  charge  are  given  by 

Vj  5=  9i^ll  +  ^2^12  *  '  '  '  ' 

^2  =  9/21  +  92^22  +  '  •  '  +(l„P2r,  > 

# 

+  92^«2  +  •  •  •  , 

where  is  the  charge  and  V.  the  voltage  on  the  ixh  wire.  In  matrix  notation,  the  above  set 

of  equations  is 

(VI  =  [PM  . 

Here  [q]  is  a  column  vector. 
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B.  Simplified  Equations 

The  above  equations  can  be  solved  for  the  general  case  of  different  voltages,  wire 
diameters,  and  locations.  This  may  be  required  in  some  cases,  for  example,  to  obtain 
charge  densities  on  multiphase  power  lines  in  proximity  to  each  other  and  to  ground.  See, 
for  example,  Doyle  et  al.  (1987),  as  well  as  Miller  (1956),  and  Adams  (1955).  However, 
these  equations  can  be  greatly  simplified  for  the  case  of  a  single  bundle  composed  of 
wires  all  having  the  same  radius  (a),  at  the  same  potential  (V),  and  located  at  a  large 
distance  above  the  earth  with  respect  to  the  bundle  dimension  h  »  b).  Stating  this  in 
terms  of  equations, 

V  =  V  =  •••  =:V=V 

a,  =  Oj  =  •  •  ■  =  a„  =  a  . 

and  the  fact  that  h  »  b  leads  to  the  approximation  that 

D^j  -  2h  =  D  for  all  i  and  j,  and 

=  92  =  •  ■  •  =  9„  =  9  • 

Substituting  the  above  into  the  matrix  equations  and  using  reciprocity  (P^  =  /V, ),  leads 
to  a  single  equation  for  the  charge  density  on  the  individual  wires 

9  =  V/(P,,  +  />,„)  . 

1.5.2  Capacitance  -  Equivalent  Radius 

The  total  charge  on  the  bundle  is 
Q  =  nq  , 

and  the  total  capacitance  of  the  bundle  is 

C  =  QIV  ^  n/(/>„  ^  P,,^  •  •  •  *  PJ  . 

A.  Single  Wire 

For  a  single  wire,  the  capacitance  is  simply  given  by 

C,  =  IjieJ  ln(2h/a)  .  (1) 

The  electric  field  at  the  surface  of  the  single  isolated  wire  obtained  using  Gauss'  law  is 
given  by 

£,  =  q/(2jt€oa)  =  ni/a}[l/ln(2A/a)l  .  (2) 

Note  the  surface  electric  field  formula  contains  one  term  related  to  geometry  enclosed 
by  {  },  and  another  term  related  to  capacitance  enclosed  by  (  ]. 
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B.  General  Case  (n-Wires) 

The  general  case  of  n  wires  uniformly  spaced  on  a  circle  of  radius  b  is  illustrated  in 
figure  1.  The  spacing  between  the  centers  of  wires  1  and  j  {dy)  is 

dy  =  2b  sin(/7t/n)  . 

Substitution  into  (1)  gives 

C„  -  2n7t€Q/{\n(2h/a)  +  ln(2/i/(2/jsin(;r/«))  + 

ln(2/i/(dsin(2jr//j))  +  •  •  ■  +  ln(2A/(2dsin((n  -  l);r//i))} 

Combining  the  log  terms  in  the  denominator  gives 

Cn  -  2n:i€Q/{ln((2h)'' /[a(2b)^"~^^  sin(jt/n)  s'm(27t/n)  *  •  ■ 

•  •  •  sm{2jn/n)  •  •  •  sin((/i  -  l)jr//i)])}. 

Further  reduction  using  the  following  identity  (appendix  A), 

2&\nin/n)  2sin(2w/«)  *  *  •  2s\n{jnfn)  •  •  * 

•  •  ■  2  sin((n  -  l)w/n)  =  n  ,  gives 
C„  =  2n7i€.J{n  ln(2b/a')}  .  or 

C„  =  2;reoAln(2/i/a')}  ■  (3) 

Where  a'  =  is  the  cage  equivalent  radius  for  capacitance,  or  the  radius 

of  a  single  wire  having  the  same  capacitance  per  unit  length  as  the  cage.  Note  that  equa¬ 
tion  (3)  has  the  same  form  as  equation  (1). 

This  simple  expression  for  the  capacitance  per  unit  length  is  quite  accurate  as  will  be 
shown  in  section  2.0.  This  result  is  equivalent  to  one  given  by  Grover  (1928).  However,  he 
did  not  obtain  the  closed  form.  The  closed  form  expression  is  included  in  gradient  formu¬ 
las  given  by  Crary  (1932)  and  Chaudhari  (Miller,  1956).  The  above  expression  for  equiva¬ 
lent  radius  is  the  same  as  the  geometric  mean  radius  of  the  cage  (Grover,  1962)  and  gives 
the  correct  inductance  in  the  high-frequency  limit.  Shelkunoff  (1952)  has  derived  the 
same  expression  for  equivalent  radius  by  a  completely  different  method  and  shows  that  it 
gives  the  correct  value  for  the  transmission  line  impedance  of  a  multiple  wire  cage. 

C.  Cylindrical  Geometry 

If  the  cage  is  located  at  the  center  of  a  cylinder  of  radius  R  instead  of  a  height  h  above 
the  ground,  the  formula  for  capacitance  uses  the  same  cage  equivalent  radius  a'  but  will 
replace  2h  by  R.  In  section  1.5.3,  capacitance  is  always  part  of  the  formula  used  for 
calculating  the  maximum  surface  electric  field  on  the  cages.  For  cylindrical  geometry,  use 
the  same  formulas  but  replace  2h  by  R  in  the  capacitance  term. 
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D.  Capacitance  In  The  Presence  Of  Other  Wires 
1.  Two-Wire  Case 

Often,  it  is  desirable  to  design  a  multiwire  cage  for  each  of  two  or  more  widely  sepa¬ 
rated  wires  carrying  a  high  voltage.  The  presence  of  the  other  wires  modifies  the  capaci¬ 
tance  of  any  particular  wire.  For  certain  conditions,  the  effect  of  this  modification  can  be 
accounted  for  by  changing  the  apparent  height  of  the  wire  . 

The  equations  for  the  two-wire  case  are 
27[€oV,  =  +  92 ‘"(f] 


where  D  is  the  slant  distance  between  one  wire  and  the  image  of  the  other  wire,  and  S  is 
the  distance  between  wires.  If  each  of  the  individual  wires  were  caged,  then  the  radius  (a) 
would  be  replaced  by 

If  both  wires  have  the  same  radius  a  (or  voltage  (hence,  the  same 

charge),  then  the  equations  reduce  to 


The  capacitance  of  each  individual  wire  is  given  by 


From  the  above  equation,  it  is  clear  that  the  capacitance  is  the  same  as  that  of  a  single 
isolated  wire  at  a  different  height.  This  new  height  is  defined  to  be  the  apparent  height 
and  is  given  by 


Thus,  for  this  case,  each  wire  can  be  treated  as  a  single  isolated  wire  at  a  new  appar¬ 
ent  height.  Note  that  the  assumptions  required  are  that  both  wires  have  the  same  (1) 
radius  a,  (2)  height  above  ground,  and  (3)  voltage.  For  the  case  of  caged  wires,  a  is 
replaced  by  which  must  be  the  same  for  both  wires. 

2.  Infinite  Wire  Grid 

The  case  of  an  infinite  number  of  parallel  wires  separated  by  a  distance  5  and  located 
at  the  same  height  h  above  the  ground  can  be  treated  in  the  same  manner.  The  capaci¬ 
tance  per  unit  length  of  one  wire  in  the  grid  is  derived  in  appendix  B  and  is  given  below; 
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where  X  is  defined  as 


This  equation  is  similar  to  that  for  the  two-wire  case.  From  this  equation,  it  can  be 
seen  that  the  concept  of  an  apparent  height  is  valid  for  an  infinite  wire  grid.  For  this  case, 
it  is  clear  that  is  given  by  the  following, 


Note  that  for  cases  like  those  above,  where  shielding  is  provided  by  other  wires  at  the 
same  potential,  the  apparent  height  is  considerably  greater  than  the  actual  height. 

1.5.3  Surface  Electric  Field 
A.  Simple  Geometric  Theory 

The  electric  field  is  calculated  by  superimposing  the  fields  due  to  the  charges  on  each 
wire.  The  field  has  a  maximum  on  the  surface  of  each  wire  at  the  point  opposite  the  cage 
center.  Figure  2  illustrates  the  method  of  calculation  of  the  maximum  field  on  one  wire  by 
vector  addition  of  the  field  due  to  the  charge  on  the  wire  itself,  plus  the  radial  component 
of  the  fields  due  to  the  charge  on  the  other  wires.  The  field  due  to  the  charge  on  the  wire 
itself  is 

E  =  ql{23t€^)  , 

where  q  is  the  charge  on  each  individual  wire,  given  by 
q  =  C„  V/n  . 

To  this  must  be  added  the  fields,  due  to  the  charges  on  the  other  wires.  The  magnitude 
of  these  fields  can  be  approximated  by  calculating  them  the  point  corresponding  to  the 
center  of  the  wire  in  question,  as  if  the  wire  were  not  there.  The  effect  of  the  wire  (cylin¬ 
der)  is  approximated  by  assuming  that  the  field  is  uniform  and,  therefore,  doubles  at  the 
point  where  the  field  is  perpendicular  to  the  surface  (Morse  &  Feshback,  1953).  Crary 
(1932)  seems  to  be  the  first  to  have  suggested  this  approximation.  The  doubling  of  the 
field  on  the  surface  at  that  point  can  also  be  derived  using  image  theory  (section  2, 
appendix  C). 
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Figure  2.  Field  calculation. 


From  symmetry  (figure  2),  it  is  seen  that,  at  the  location  of  maximum  field  on  the 
surface  of  a  wire,  the  field  due  to  the  charge  on  the  other  wires  will  be  normal  to  the 
surface.  If  the  field  were  uniform,  then  the  effect  of  the  cylinder  on  this  field  can  be 
accounted  for  by  multiplying  by  a  factor  of  2.  The  uniform  field  approximation  is  accu¬ 
rate  when  the  wires  are  well  separated,  but  becomes  less  accurate  when  the  wires  are 
closer  together.  The  accuracy  of  this  approximation  will  be  examined  further  in  the  next 
section. 

From  figure  1  and  Gauss’  law,  it  can  be  seen  that  the  field  at  the  location  of  the  center 
of  wire  1  due  to  the  charge  on  the  ^  wire,  is  given  by 

E  =  q/(2jt€Qd^p  =  SinQn/n))  . 

To  get  the  component  in  the  direction  of  the  cage  radial  requires  multiplication  of  the 
above  term  by  cos(3Jy),  which  cancels  out  the  smijidn)  term  (figure  D-1,  appendix  D). 
Thus,  the  contribution  to  the  radially  directed  electric  field  at  the  location  of  the  center  of 
a  wire,  due  to  the  charge  on  each  of  the  other  wires,  is  identical  and  given  by 


E  *  qliAne^b)  . 

The  maximum  surface  electric  field  then  is  approximated  by  the  sum  of  the  field  from 
the  wire  itself  plus  twice  the  sum  of  the  contributions  from  the  other  wires,  given  by 
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Substituting  for  q  by  using  the  capacitance  formula  ,  gives 


Substituting  for  gives 

where  { }  contains  the  geometry  term,  and  [  ]  the  capacitance  term.  Equations  (3)  and  (4) 
are  called  the  simple  geometric  theory  (SGI^.  This  formula  appears  in  Chaudhari  (Miller, 
1956)  and  Smith  (1985). 

B.  Complete  Geometric  Theory 

The  simple  geometric  formula  is  based  on  estimating  the  field  in  the  vicinity  of  one 
cylinder  due  to  all  the  other  cylinders  by  calculating  the  field  at  the  center  of  the  wire  as  if 
it  were  not  there.  This  field  is  then  doubled  to  account  for  the  effect  of  the  cylinder.  This 
would  be  a  valid  approach  if  the  field  in  the  vicinity  of  the  cylinder  were  uniform.  This 
assumption  is  violated  when  the  wires  are  in  close  proximity.  From  the  discussion  of  the 
image  method  in  appendix  C,  it  is  clear  that  a  better  approximation  would  be  to  use  twice 
the  field  calculated  at  the  surface  of  the  cylinder  instead  of  at  the  center  of  the  cylinder. 

The  complete  geometric  theory  (CGT)  is  derived  by  the  same  method  as  the  simple 
geometric  theory  (SGT)  except  that  the  fields  due  to  Ae  charge  on  the  other  wires  are 
calculated  on  the  surface  of  the  wire,  not  at  the  center  of  the  wire.  The  derivation  is  given 
in  detail  in  appendix  D.  The  simple  formula  for  capacitance  (equation  3)  is  used  in  this 
derivation. 

The  resulting  equation  for  the  maximum  surface  field  is  similar  to  SGT  but  more 
complicated.  The  geometric  factor  for  the  two-wire  cage  is 

The  geometric  terms  for  three  or  more  wires  are  increasingly  more  complicated.  They 
are  derived  explicitly  in  appendix  D  for  the  three-  and  four-wire  case.  The  general  expres¬ 
sion  for  the  geometric  term  is  a  finite  summation,  which  is  also  given  in  the  appendix. 

C.  Modified  Simple  Geometric  Theory 

Note  that  for  the  two-wire  case,  the  only  difference  between  SGT  and  CGT  is  that  the 
factor  a/2  has  been  added  to  the  denominator  of  the  second  term  in  the  geometric  factor. 
This  suggests  a  modification  to  the  SGT  formula  having  a  geometric  factor  of  the  form 

Replacing  the  geometric  factor  in  equation  (4)  with  the  above  factor  results  in  the 
following  equation,  which  is  called  the  modified  simple  geometric  theory  (MSGT): 
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(5) 


E  =  ^-iir  1  1 

"  n  ^  6  +  fl/2jLln(2A/a')J 

Appendix  D  also  contains  an  alternate  derivation  of  the  MSGT  formula. 

As  will  be  shown  in  section  2.0,  the  MSGT  formula  provides  a  surprisingly  good  fit  to 
reality  for  all  cases  considered. 

Note  that  for  the  two-wire  case,  the  formulas  for  the  MSGT  and  CGT  are  the  same. 
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2.0  COMPARISON  WITH  MORE  EXACT  METHODS 


2.1  TWO-WIRE  CASE 

This  section  compares  the  results  of  the  SGT  and  MSGT  to  several  other  methods  of 
calculation  for  two- wire  cages. 

Substitution  of  n  =  2  in  equations  (3),  (4),  and  (5)  gives 
Cj  -  2jrc„[l/ln  ■ 

l(i  ^  it'"  (p^)] 

•  Hi  ^  i^i«}['^'"  ((^)] 


2.1.1  Charge  Offset  Method 

Appendix  E  contains  a  different  method  for  the  derivation  of  capacitance  and  surface 
gradient  for  the  two-wire  case.  This  derivation  uses  a  single  line  charge  near,  but  offset 
from,  the  center  of  each  wire.  For  the  case  of  two  wires  with  equal  charge  of  opposite 
sign,  this  is  an  exact  solution.  However,  for  the  multi  wire  cage  geometry  that  we  are 
using,  where  the  wires  are  at  the  same  potential,  this  method  is  approximate. 

For  this  case,  the  offset  is  adjusted  to  give  equal  potential  at  the  wire  surface,  at  two 
points  located  on  opposite  sides  of  the  wire.  It  is  interesting  to  note  that  the  offset  cancels 
out  implicitly,  and  the  resulting  formulas  for  capacitance  and  surface  gradient  are  exactly 
the  same  as  the  SGT. 

2.1.2  Image  Theory 

The  electric  field  in  the  vicinity  of  a  bundle  of  charged  conductors  can  be  calculated 
by  using  the  method  of  multiple  image  charges  (EPRI,  1987).  A  charge  is  assumed  on  one 
cylinder  and  imaged  on  the  other  cylinder.  The  image  charge  is  then  imaged  back  into  the 
first  cylinder,  etc.  This  leads  to  an  infinite  series  solution,  which  is,  in  practice,  truncated. 
Care  must  be  taken  to  construct  a  series  that  converges  rapidly.  The  secret  is  to  develop  a 
series  of  images  using  pairs  of  charges  forming  dipoles,  that  have  no  net  charge,  and  ever 
decreasing  dipole  moment  (appendix  C).  The  series  converges  to  within  0.1%  with  40 
terms  when  the  two  wires  are  touching,  and  much  faster  as  the  two-wires  are  spaced 
farther  apart.  Computer  calculations  were  done  for  various  values  of  2/i/a,  using  100 
terms  of  the  series.  The  results  of  these  calculations  are  given  in  section  2.1.5  and  are 
compared  to  the  SGT  and  MSGT. 
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For  the  general  n-mrt  cage,  the  image  formulation  becomes  increasingly  cumbersome 
because  each  charge  gives  rise  to  an  image  on  each  of  the  other  cylinders,  which  gives 
rise  to  a  second  image  on  each  of  the  other  cylinders,  etc.  The  location  of  the  images 
becomes  complicated  in  that  they  are  not  all  located  on  a  single  radius  from  the  wire 
center  as  they  are  for  the  two-wire  case.  For  this  reason,  the  image  theory  was  not  pur¬ 
sued  for  n  larger  than  2.  If  there  was  no  other  alternative,  image  theory  would  be  tedious 
but  useful.  However,  the  method-of-moments  approach  described  below  has  advanced  to 
the  point  where  it  is  easy  to  apply. 

2.1.3  Conformal  Transformation  Theory 

The  field  distortion  from  nearby  wires  for  the  two-wire  case  can  also  be  solved  by 
using  the  theory  of  conformal  transformations  (Poritsky  in  Clark,  1932).  This  also  leads  to 
an  infinite  series  that  must  be  truncated  for  calculations.  This  technique  was  used  in  the 
reference  to  evaluate  the  assumption  that  the  field  distortion  caused  by  the  proximity  of 
the  cylinders  is  the  same  as  if  the  field  in  the  vicinity  of  the  cylinders  is  uniform.  The 
conclusion  was  that  the  assumption  gives  results  that  are  accurate  when  the  wires  are 
separated  by  more  than  5  wire  diameters.  Conformal  transformation  and  imag  -theory 
results  are  numerically  the  same. 

2.1.4  Method  Of  Moments 

The  method  of  moments  (MOM)  is  a  computer  technique  for  solving  a  general  set  of 
linear  equations  (Harrington,  1968).  In  this  case,  the  surface  of  the  wires  is  divided  into 
many  segments,  and  the  charge  on  each  segment  determined  to  match  the  boundary  con¬ 
ditions.  If  enough  segments  are  used,  the  solution  can  be  considered  to  be  exact  as  it 
accounts  for  field  distortion  due  to  the  proximity  of  the  wires  to  each  other  and  to  ground. 

A  commercially  available  computer  code  was  used  to  calculate  capacitance  and  elec¬ 
tric  field  for  the  two-wire  case  (H£CTRO).2  The  results  for  the  two-wire  case  are  com¬ 
pared  to  the  other  methods  of  calculation  in  the  following  section.  As  an  initial  test,  the 
MOM  code  was  exercised  for  several  cases  of  a  single  wire  above  ground  and  found  to 
give  numerically  accurate  results  when  the  number  of  segments  per  wire  was  16  or  more. 

2.1.5  Comparison  Of  Results  (Two- Wire  Case) 

Figure  3a  shows  a  comparison  of  calculated  values  of  capacitance  (Q)  for  the  two- 
wire  cage  with  2h/a  =  5000,  The  results  have  been  normalized  by  division  by  the  capaci¬ 
tance  of  a  single  wire.  Calculations  were  done  using  MOM,  image  theory,  and  SGT. 
Figure  3b  gives  the  percentage  error  for  the  MOM  and  SGT  by  using  the  MOM  calcula¬ 
tions  as  the  reference.  Note  that  the  MOM  results  are  numerically  the  same  as  the  image 
theory  except  when  the  wires  are  touching.  The  SGT  results  are  within  1.5%  when  the 
wires  are  touching  and  better  than  0,5%  when  bla  >  2. 


*  2.  ELECTRO-a  computer  program  available  from  Integrated  Engineering  Software,  Winnipeg,  Manitoba,  Canada. 
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CAPACITANCE  PERCENTAGE  ERROR  NORMAUZED  CAPACITANCE 
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Figure  3a.  Normalized  capacitance— two-wire  cage 
{2hla  »  5000). 
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Figure  3b.  Capacitance  percentage  error— two-wire  cage 
{2hla  =  5000)  (MOM). 
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Figure  4a  gives  the  maximum  surface  electric  field  for  the  same  case.  The  calculations 
were  done  by  using  the  same  methods  as  above  as  well  as  CGT.  Note  that,  for  the  two- 
wire  case,  the  CGT  formula  is  the  same  as  the  MSGT  formula.  Again,  the  results  have 
been  normalized  to  the  value  for  a  single  wire.  Figure  4b  gives  the  percentage  error  for 
E-field,  with  MOM  as  the  reference.  Again,  note  that  the  MOM  results  are  numerically 
identical  to  the  image-theory  results.  Note  that  SGT  does  not  give  good  results  for  small 
values  of  bla.  However,  CGT  gives  excellent  results  for  all  values  of  b/a  except  when  the 
wires  are  touching. 

From  the  figures,  it  is  concluded  that  MOM  gives  numerically  exact  results  except 
when  the  wires  touch.  The  MOM  results  will  be  used  as  the  reference  for  cases  involving 
three  or  more  wires. 

Another  conclusion  is  that  the  SGT  capacitance  results  are  very  accurate  for  all  values 
of  b/a.  The  maximum  error,  less  than  2%,  occurs  when  the  wires  touch. 

For  the  two-wire  case,  the  CGT  gives  excellent  predictions  for  surface  gradient. 

Figures  5a  and  b,  and  6a  and  b  are  similar  to  the  above  but  with  2h/a  =  300.  The  SGT 
capacitance  values  are  again  very  accurate  for  small  values  of  b/a.  However,  in  this  case, 
they  diverge  at  larger  values.  This  is  because  the  approximation  that  =  2h  is  poor. 
Because  of  this  assumption,  the  SGT  capacitance  formula  is  valid  only  when  h  >  5b. 

2.2  GENERAL  CASE,  w-WIRE  CAGE 


2.2.1  MOM  Results  For  3,  4,  5,  6,  and  12  Wires 

The  MOM  has  been  validated  in  the  above  section  for  the  geometry  involved  in  multi¬ 
wire  cages.  The  validation  used  the  case  of  a  two-wire  cage  and  validated  MOM  against 
exact  methods  of  calculation.  When  there  are  three  or  more  wires,  the  exact  image-theory 
calculations  become  very  complicated.  The  only  practical  “exact”  method  of  calculation  is 
to  use  the  method  of  moments. 

To  examine  the  accuracy  of  the  simple  formulas  (SGT,  MSGT,  and  CGT)  for  n  greater 
than  2,  a  data  set  was  developed  by  using  MOM.  The  cases  run  were  «  =  3,  4,  5,  6,  and 
12,  with  2h/a  *=  5000.  The  results  for  capacitance  and  surface  electric  field  are  given  in 
figures  7  and  8.  The  two-wire  case  has  been  included  in  the  figures  for  completeness. 

Note  that,  due  to  the  geometry  and  the  “exact"  nature  of  the  solution,  the  charge  on 
each  wire  is  not  the  same.  The  charge  on  the  wires  nearest  to  the  image  is  slightly  greater 
than  the  charge  on  wires  farther  away  from  the  image.  Thus,  the  electric  field  and  capaci¬ 
tance  is  slightly  different  depending  upon  which  wire  is  examined.  For  this  reason,  the 
plots  in  the  figures  include  a  maximum,  a  minimum,  and  the  average  value  for  both 
capacitance  and  surface  electric  field. 
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SURFACE  E  FIELD.  PERCENTAGE  ERROR  MAXIMUM  SURFACE  E  FIELD  (NORMAUZED) 
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Figure  4a.  Normalized  surface  E-field— two-wire  cage 
{Ihla  X  5000). 
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Figure  4b.  E-field  percentage  error— two-wire  cage 
(2/i/a  X  5000)  (MOM). 
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CAPACfTANCE  PERCENTAGE  ERROR  NORMAUZED  CAPACfTANCE 


b/a 

Figure  5a.  Normalized  capacitance— two-wire  cage 
{2hla  =  300). 


Figure  5b.  Capacitance  percentage  error— two-wire  cage 
(2/j/fl  *=  300)  (MOM). 
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SURFACE  E  FIELD,  PERCENTAGE  ERROR  MAXIMUM  SURFACE  E  FIELD  (NORMAUZEO) 
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Figure  6a.  Normalized  surface  E-field— two-wire  cage 
(2/i/a  =  300). 
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Figure  6b.  E-field  percentage  error— two-wire  cage 
{Ihla  =  300)  (MOM). 
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2.2.2  Capacitance  Comparison 

The  simple  geometric  theory  for  capacitance,  as  given  in  equation  (3),  has  been  com¬ 
pared  with  the  value  of  capacitance  calculated  by  the  MOM.  The  results  are  presented  in 
figures  9  through  13.  The  results  are  presented  if  the  form  of  plots  of  the  ratio  of  the  SGT 
calculated  value  to  the  “actual”  value  are  as  calculated  by  using  the  MOM  results.  The 
comparison  has  been  made  to  the  maximum,  minimum,  and  average  values  calculated  by 
MOM. 

These  figures  show  that  the  value  of  capacitance,  by  using  SGT,  is  always  within  2.5% 
of  the  actual  average  value.  For  values  of  bla  greater  than  3,  the  SGT  gives  values  that  are 
within  1%  of  the  actual  value. 

The  separation  between  maximum  and  minimum  values  increases  with  separation  but 
in  general  is  less  than  1%  for  all  values  of  bla  considered. 

The  conclusion  is  that  the  SGT  capacitance  formula  is  an  excellent  approximation! 

2.2.3  Surface-Electric-Field  Comparison 

Three  formulas  for  maximum  surface  electric  field  are  to  be  examined,  SGT,  MSGT, 
and  CGT.  The  comparison  has  been  done  by  taking  the  ratio  of  the  result  calculated  by 
the  formulas,  and  the  actual  value  as  calculated  by  the  MOM.  The  results  are  presented  in 
figures  14  through  18.  For  the  MSGT  formula,  the  comparison  has  been  made  to  the 
maximum,  minimum,  and  average  values  calculated  by  MOM.  The  SGT  and  CGT  have 
been  compared  only  with  the  average  value  calculated  by  the  MOM. 

The  SGT  formula  has  large  errors  for  small  values  of  b/a,  decreasing  as  bla  becomes 
larger.  For  example,  for  the  three-wire  cage  with  bla  of  1.5,  the  SGT  error  is  in  excess  of 
30%.  Because  of  the  large  error  for  small  spacings,  the  SGT  formula  gives  an  erroneous 
value  for  the  separation  having  the  minimum  surface  electric  field. 

CGT  results  have  been  calculated  for  n  =  3  and  4.  The  error  for  the  CGT  method  is 
significant  at  small  spacings  but  less  than  for  the  SGT  method. 

The  MSGT  formula  is  best  for  all  spacings.  It  has  less  than  a  6%  error  at  very  small 
values  of  bla,  and  less  than  2%  error  for  values  of  bla  of  three  or  greater. 

Hie  MSGT  formula  gives  the  most  accurate  results  of  the  simple  formulas  tried.  It 
also  gives  the  correct  value  of  separation  that  has  minimum  surface  electric  field.  This 
formula  should  be  adopted  as  the  standard  formula  for  cage  calculations. 
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C^(sgt)/C^(,MOM)  C3(MSGT)fC3(MO\f) 


Figure  9.  Capacitance  comparison  (MSGT/MOM)— three-wire 
cage  (2h/a  =  5000). 
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Figure  10.  Capacitance  comparison  (SGT/MOM)— four-wire 
cage  (2h/a  =  5000). 
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Figure  11.  Capacitance  comparison  (SGT/MOM)— five-wire 
cage  {Ihia  =  5000). 
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Figure  12.  Capacitance  comparison  (SGT/MOM)— six-wire 
cage  {Ihia  =  5000). 


24 


Figure  13.  Capacitance  comparison  (SGT/MOM)— twelve-wire  cage 
{Ihia  =  5000). 


Figure  14.  Normalized  surface  E-field  comparison  ‘o  MOM 
three-wire  cage  (2/i/a  =  5000). 
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Figure  15.  Maximum  surface  E-field  comparison  to 
MOM— four-wire  cage  (2h/a  =  5000). 
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Figure  16.  Maximum  surface  E-field  comparison  to 
MOM— five- wire  cage  (2h/a  =  5000). 
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Figure  17.  Maximum  surface  E-field  comparison  to 
MOM— six-wire  cage  (2A/a  «  5000). 
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Figure  18.  Maximum  surface  E-field  comparison  to 
MOM— twelve-wire  cage  (2A/a  «  5000). 
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3.0  DESIGN 


3.1  DESIGN  DATA 

Design  curves  have  been  developed  by  using  the  MSGT  formula. 

Figure  8  shows  that,  at  first,  the  surface  gradient  decreases  rapidly  with  increasing  b, 
then,  more  slowly,  goes  through  a  broad  minimum  and  then  slowly  increases. 

The  reason  for  the  minimum  is  as  follows:  The  surface  gradient  is  the  product  of  a 
capacitance  factor  and  a  geometry  factor.  The  capacitance  factor  increases  with  separa¬ 
tion  while  the  geometry  factor  decreases  with  separation.  When  wires  are  less  than  one 
wire  radius  apart,  the  product  is  somewhat  less  than  the  value  for  a  single  isolated  wire. 
At  large  separations,  both  factors  and  the  product  are  asymptotic  to  the  value  for  a  single 
isolated  wire.  For  close  separations,  the  geometry  factor  falls  off  much  faster  than  the 
capacitance  increases  and,  thus,  the  product  of  the  two  factors  has  a  minimum  value. 

The  capacitance  factor  is  a  logarithmic  function  of  height.  Thus,  the  rate  increase  of 
the  capacitance  is  reduced  as  height  increases.  The  rate  decrease  for  the  geometric  factor 
is  independent  of  height.  Thus,  the  location  of  the  minimum  is  a  function  of  height.  The 
minimum  occurs  at  wider  spacings  as  height  increases. 

A  computer  program  based  on  MSGT  was  developed  to  find  the  value  of  b/a,  for 
which  the  surface  gradient  is  minimum.  This  program  is  included  in  appendix  F.  The 
location  of  the  minimum  is  called  b^/a  and  is  plotted  versus  2h/a  in  figure  19.  From 
figure  19,  bj^/a  is  seen  to  increase  approximately  logarithmically  with  height,  as  well  as 
linearly  with  n. 

The  wire-to-wire  spacing  at  minimum  gradient  is  called  5^,  and  is  plotted  in  figure 
20.  From  this  figure,  it  is  seen  that  is  nearly  the  same  for  n  =  3  to  6.  also  varies 
logarithmically  with  height.  Curiously,  the  two-wire  case  has  a  different  slope  than  all  the 
other  cases.  As  n  increases,  the  slope  of  5^^  remains  nearly  the  same  as  cases  n  =  3,  but 
the  curve  moves  down  slightly  with  n.  This  is  illustrated  in  figure  20  by  the  curve  for  n  = 
12. 

Table  1  gives  values  of  in  terms  of  wire  diameters  for  practical  heights.  For 
1 -inch-diameter  wires,  the  corresponding  heights  in  the  table  are  20.8,  208,  and  2083  feet. 
Table  1  shows  that  for  minimum  gradient,  wire-to-wire  spacings  need  not  be  more  than  8- 
to  15-wire  diameters  for  the  two-wire  case  and  slightly  more  when  n  is  greater  than  2. 


Table  1.  S^I2a  for  practical  values  of  2h/a. 


n 

2hla 

2h/a 

2hla 

1000 

10,000 

100,000 

2 

8.86 

13.11 

17.45 

3-6 

9.75 

15.5 

21.5 

3.1.1  Equivalent  Radius  for  Electric  Field 

The  concept  of  equivalent  radius  for  electric  field  of  a  cage  is  a  useful  design 
concept  that  was  introduced  by  Miller  (1956).  This  concept  is  defined  as  the  radius  of  a 
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single  wire  that  would  have  the  same  maximum  surface  gradient  as  the  cage  when  located 
at  the  same  position  as  the  cage  axis. 

It  is  defined  by  the  equation 


Note  that  this  equation  has  the  same  form  as  equation  (2). 

Using  normalized  parameters  gives  the  following  equation  for  , 


£n, 

a 

1 

3* 

^1 

^«qE 

in 

VO  “.tcJ 

This  transcendental  equation  can  be  solved  by  iterative  search  or  plotting.  An  iterative 
search  is  included  in  the  computer  program  that  finds  to  give  the  corresponding 
value  of  A  sample  of  the  output  of  this  program  is  included  in  appendix  F,  The 

values  of  for  b  =  b^.„  are  useful  for  design  and  are  plotted  in  figure  21. 

It  is  interesting  that  the  abscissa  of  figure  21  can  be  converted  from  2h/a  to  by 

dividing  the  abscissa  values  by  the  value  of  the  function  This  gives  a ^  as  a 

function  of  2hla^g^  for  cages  with  radius  b  *  b^.  A  set  of  tfiese  curves  for  values  of  n 
from  2  through  6  and  ^  «  12  are  given  in  figure  22.  These  curves  are  even  more  useful  for 
design. 

Finally,  it  was  observed  that  the  curves  of  figure  19  could  be  normalized  to  by 
dividing  ^e  values  of  both  the  ordinate  and  abscissa  by  the  value  of  the  function  a^g^/a  at 
each  point.  The  resulting  set  of  curves,  given  in  figure  23,  completes  the  design  data. 
Once  a  value  for  has  been  determined,  the  curves  of  figures  22  and  23  can  be  used 
to  determine  the  optimum  cage  dimensions  meeting  the  surface  electric  field  requirement. 

Note  that  the  curves  of  figure  23  are  nearly  straight  lines  on  the  linear-log  scale.  All 
the  curves  have  nearly  the  same  slope.  These  curves  can  be  approximated  by  the  follow¬ 
ing  simple  formula: 

^min  =  2.36  X  Log,ola,,£/a]  -  /„ , 
where  /„  is  the  intercept  value  and  is  given  in  table  2. 

Table  2.  Intercept  values  for  b  ■  equation. 


n 

K 

2 

0.8939 

3 

0.7184 

4 

0.6095 

5 

0.5216 

6 

0.4684 

12 

0.2982 

30 
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Fif  ure  22.  la  versus  2hla^^^  for  various  values  of  n. 
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10  100  1000  10000  100000 
Figure  23.  versus  2hla^^^  for  various  values  of  n. 


3.2  DESIGN  PROCEDURES 

A  major  design  consideration  for  cages  is  to  minimize  the  structural  impact  on  the  rest 
of  the  structure.  Increasing  the  spacing  of  the  cage  increases  the  projected  wind  area,  and 
the  spreader  hardware  weight.  Cost  considerations  and  wire  availability  often  dictate  the 
size  of  wires  available  to  the  designer.  For  these  reasons,  the  first  design  procedure  starts 
with  a  given  diameter  of  wire  to  be  used  in  the  cage  and  determines  the  minimum  spacing 
that  will  have  a  satisfactory  gradient. 

3.2.1  Minimum  Spacing 

For  this  case,  the  following  parameters  are  given: 

height  h 
wire  radius  a 
design  voltage  V 
design  gradient  £o 

The  design  procedure  consists  of  finding  the  value  of  the  cage  radius  b  that  has  a 
gradient  equal  to  £o-  The  approach  is  to  find  E[^Ey  using  the  single  wire  formula,  and 
then  plot  equation  (5)  for  several  values  of  n.  The  design  value  of  b  is  selected  as  the 
value  where  £„/£i  =  Safety  factors  are  included  in  the  selection  of  Eq.  To  mini¬ 

mize  the  structural  loading  effects,  the  minimum  value  of  n  is  chosen  that  satisfies  the 
gradient  criteria. 
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A.  First  Design  Example  (Given  Wire  Size) 

As  an  example,  select  1 -inch-diameter  wires,  400  feet  above  the  ground,  operated  at 
200  kV  with  a  design  gradient  of  1  kV/mm.  The  value  of  2hla  -  19,200.  The  design 
parameter  is 

Eo/fi  =  0.6263  =  /£i  . 

Figure  24  has  plots  of  E^IE-\  for  several  values  of  n. 

Note  that  the  «  =  1  curve  is  a  special  case  for  the  single  wire.  It  is  the  ratio  of  £  for  a 
wire  with  radius  b,  to  £  for  a  1 -inch-diameter  wire.  This  curve  can  be  used  to  find 
which  is  the  ordinate  value  where  £„/£i  =  0.6263. 

From  figure  24,  the  following  conclusions  can  be  drawn: 

(1)  The  n  =  1  case  gives  =  1.68. 

(2)  For  1 -inch-diameter  wires,  there  is  no  spacing  that  will  give  an  acceptable 
gradient  for  n  =  2. 

(3)  For  n  =  3,  the  design  criteria  will  be  met  when  bla  =  2.4. 

(4)  For  n  *  4,  the  cage  radius  is  only  slightly  more  than  that  of  <.  single  wire  meeting 
the  criteria.  For  higher  values  of  n,  the  minimum  radius  of  the  cage  will  be 
larger  than  the  radius  of  a  single  wire  that  meets  the  criteria. 

Therefore,  the  best  cage  meeting  the  criteria  is  a  three-wire  cage  with  a  cage  radius  of 
1.2  inches  using  1-inch  wires. 


Figure  24.  First  design  example,  EJE\  for  various  values  of 
n  {Ihia  =  19,200). 
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3.2.2  Minimum  Wire  Cross  Section 

For  the  second  design  procedure,  the  wire  diameter  is  not  specified  but  will  be  deter¬ 
mined.  The  design  parameters  are  the  same  as  above,  but  wire  radius  is  not  specified. 
This  design  procedure  results  in  a  cage  with  the  minimum  cross  sectional  area  that  meets 
the  gradient  design  criteria.  It  starts  with  determination  of  the  required  radius  for  the 
single  wire  that  will  meet  the  design  criteria.  This  is  the  required 

Once  is  known,  calculate  the  parameter  and  then  is  determined  by 

figure  22  for  each  value  of  n  chosen.  The  required  value  of  wire  radius  for  the  cage  is 
determined  by 

a  = 


Next,  determine  b^/a  from  figure  23  for  each  value  of  n. 

The  above  procedure  defines  the  cage  parameters  b  and  a.  As  a  check,  use  these 
parameters  in  equation  (5)  to  be  sure  the  design  criteria  is  met. 

A.  First  Design  Example  (Minimum  Wire  Cross  Section) 

This  example  uses  the  same  design  parameters  as  in  section  3.2.1. 

Design  Parameters 
Height  h  *  400  feet 
Voltage  V  »  200  kV 
£o  =  IkV/mm 

From  the  above  example,  we  already  know  that  the  required  radius  of  a  single  wire  is 
=  1-68  X  0.5  inch  =  0.84  inch. 

This  gives  a  value  of  2hla^^^  =  11,429.  From  the  curves  in  figures  22  and  23,  deter¬ 
mine  and  b^la.  The  design  procedure  is  demonstrated  by  table  3.  The  results  are 

progressively  obtained  moving  from  left  to  right  in  the  table.  The  notation  at  the  bottom 
of  each  column  indicates  how  the  results  in  that  column  were  derived.  The  following 
conclusions  follow  from  table  3. 

(1)  The  minimum  required  wire  diameter  for  a  two-wire  cage  is  2  x  0.515  =  1.03 
inches.  This  is  consistent  with  the  result  of  the  first  part  of  this  example  that 
indicated  that  there  was  no  two-wire  solution  for  1-inch-diameter  wires. 

(2)  The  minimum  cross-section  cage  radius  is  about  the  same  for  all  values  of  n. 

(3)  An  optimum  solution  can  be  found  for  any  value  of  n.  The  procedure  outlined 
finds  these  solutions.  The  £„  column  is  the  calculated  value  of  maximum  gradi¬ 
ent  using  equation  (5).  Note  that  the  design  criteria  is  met  exactly! 

(4)  The  area  column  contains  the  cross-sectional  area  of  the  cage  wires.  The  weight 
and  strength  of  the  cage  will  be  proportional  to  this  value.  The  total  cross- 
sectional  area  of  the  cage  wires  decreases  significantly  as  n  increases. 
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Consequently,  the  weight  of  the  cage  will  be  less  for  larger  values  of  n.  The 
weight  will  not  decrease  as  much  as  the  cross-sectional  area  because  of  the 
added  weight  of  the  spreader-hardware.  The  cross-sectional  area  decreases  by  a 
factor  of  4.12  from  1  wire  to  12  wires. 

(5)  The  windage  column  contains  the  projected  profile  for  windage.  The  profile  is 
calculated  by  assuming  each  wire  acts  independently  and  uses  a  factor  of  1.2  to 
account  for  the  cylindrical  shape  of  the  wires.^  TTiis  windage  profile  increases 
with  n,  indicating  that  wind  loading  will  be  higher  if  more  wires  are  used.  The 
windage  increases  48%  from  1  wire  to  12  wires. 


Table  3.  First  design  example-minimum  wire  cross-section  wire  diameter  and  cage 
radius  for  various  values  of  n,  200  kV,  h  =  400  feet.a^  ^  =  0.84inch,  2hla^  ^  *  11, 429, 
£o  =  1  kV/mm.  ^ 


n 

^van/^tqE 

a  inch 

d  inch 

b  inch 

^tqC 

Area 

inch2 

Windage 

inch 

1 

1 

0.84 

1.680 

1.00 

2.217 

2.016 

2 

1.63 

S.15 

0.515 

1.031 

7.350 

2.752355 

1.00 

1.669 

2.474 

3 

2.27 

8.93 

0.370 

0.740 

7.501 

3.96776 

1.00 

1.291 

2.664 

4 

2.932 

9 

0.286 

0.573 

7.560 

4.71721 

1.00 

1.031 

2.750 

5 

3.559 

9.063 

0.236 

0.472 

7.613 

5.243548 

1.00 

0.875 

2.832 

6 

4.2 

9.125 

0.200 

0.400 

7.665 

5.627174 

1.00 

0.754 

2.880 

12 

8.1 

9.25 

0.104 

0.207 

7.770 

6.670124 

1.00 

0.405 

2.987 

Fig  22 

Fig  23 

0.84/CO12 

2*a 

0.84/CO13 

Eqn  5 

n*Pi*a2 

1.2*n'a 

3.3  SECOND  DESIGN  EXAMPLE 

This  example  corresponds  to  the  end  of  span  E-5  in  the  east  array  at  VLF  Jim  Creek. 
The  end  of  this  span  is  quite  far  from  the  other  spans  and  is  treated  as  if  it  were  a  single 
isolated  wire. 

The  existing  span  is  a  copper  weld  wire  1.01  inches  in  diameter.  The  end  of  the  span 
is  600  feet  away  from  ground.  Given  voltages  of  200  kV  and  330  kV  and  a  design  gradient 
of  1  kV/mm  design  an  adequate  cage  for  this  wire. 

3.3.1  Given  Wire  Size 

First,  calculate  the  gradient  on  the  wire  for  the  two  design  voltages.  The  simple  single¬ 
wire  formula  can  be  used  because  the  end  effect  is  assumed  to  be  mitigated  by  the  corona 
rings  on  the  insulators. 

a  s  1.01  inches 
h  =  600  feet 
Eq  *  1.00  kV/mm 

*  3.  Personal  communication  with  Bob  Prince,  NAVFAC  Chief  Engineer  for  Towers. 
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Ihia  =  28,515 

V  =  200  kV  £i  =  1.52  kV/mm  EqIE^  «  0.66 

V  =  330  kV  £i  =  2.50  kV/mm  £o/£i  =  0.40 

Figure  25  has  the  corresponding  plot  of  equation  (5).  The  value  of  and  bta  versus 
n  are  obtained  from  the  figure  as  the  points  where  £„  is  equal  to  0.66  and  0.40  respec¬ 
tively.  The  results  are  given  in  tables  4a  and  4b.  The  positions  with  an  X  indicate  either 
no  solution  or  a  single  wire  will  have  the  same  or  a  smaller  total  diameter  as  the  cage. 

The  following  conclusions  can  be  drawn  from  tables  4a  and  4b. 

(1)  For  a  fixed  wire  size,  both  weight  and  windage  increase  with  n.  Therefore,  the 
cage  with  the  smallest  value  of  n  that  satisfies  the  design  criteria  is  selected. 

(2)  For  200  kV,  a  two-wire  cage  with  5-1/2  (2  x  2.75)  -inch  separation  would  be 
adequate. 

(3)  For  330  kV,  a  four-wire  cage  would  be  required. 


Table  4a.  Second  design  example— given  wire  diameter  and  cage  radius  for  various 
values  of  «,  a  =  0.505,  2h/a  =  28,515  200  kV,  =  0.803  inch. 


n 

bja 

a  inch 

d  inch 

b  inch 

Area 

Windage 

1 

1.590 

0.803 

1.606 

1.00 

2.025 

1.927 

2 

1.664203 

5.430 

0.505 

1.010 

2.74215 

1.00 

1.602 

2.424 

3 

1.125121 

1.920 

0.505 

1.010 

0.9696 

1.00 

2.404 

3.636 

4 

0.987295 

1.540 

0.505 

1.010 

0.7777 

1.00 

3.205 

4.848 

5 

X 

6 

X 

Table  4b.  Second  design  example- 

-330  kV,  a 

=  1 .49  inches. 

n 

«e,c/a 

bja 

a  inch 

d  inch 

b  inch 

Area 

Windage 

1 

2.790 

1.409 

2.818 

1.00 

6.237 

3.381 

2 

X 

3 

X 

4 

4.198488 

10.610 

0.505 

1.010 

5.35805 

1.00 

3.205 

4.848 

5 

2.713159 

5.470 

0.505 

1.010 

2.76235 

1.00 

4.006 

6.060 

6 

2.304374 

4.320 

0.505 

1.010 

2.1816 

1.00 

4.807 

7.272 

3.3.2  Minimum  Wire  Cross  Section 

By  using  the  second  design  procedure  to  find  the  cage  with  minimum  wire  cross  sec¬ 
tion,  the  above  design  constraints  will  be  satisfied  with  the  design  curves  in  figures  22  and 
23.  The  results  are  given  in  table  5  for  the  200-kV  case. 
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The  following  conclusions  follow  from  tables  4a  and  5: 

(1)  Cages  designed  for  minimum  wire  cross  section  have  less  area  and  windage  for 
the  same  value  of  n  than  cages  designed  by  the  fixed  wire  size  procedure. 

(2)  The  two-wire-minimum  cross-section  cage  satisfying  the  design  criteria  uses 
wires  with  0.98-inch  diameter  and  a  separation  of  14.8  inches. 

(3)  All  the  minimum  cross-section  cages  have  a  cage  radius  at  about  7.5  inches. 


Table  5.  Second  design  example— minimum  wire  cross-section  wire  diameter  and  cage 
radius  for  various  values  of  n,  200  kV,  h  =  600  feet.a^^^.  =  0.803  inch,  2h/a^^^  =  17, 933, 
£o  =  1  kV/mm. 


n 

a  inch 

d  inch 

b  inch 

En 

Area 

inch2 

Windage 

inch 

1 

1.000 

0.803 

1.606 

1.00 

2.026 

1.927 

2 

1.640 

9.200 

0.490 

0.979 

7.388 

2.689692 

1.00 

1.506 

2.350 

3 

2.304 

9.297 

0.349 

0.697 

7.465 

3.876951 

1.00 

1.145 

2.509 

4 

2.980 

9.477 

0.269 

0.539 

7.610 

4.668524 

1.00 

0.912 

2.587 

5 

3.610 

9.510 

0.222 

0.445 

7.637 

5.194605 

1.00 

0.777 

2.669 

6 

4.290 

9.557 

0.187 

0.374 

7.674 

5.570992 

1.00 

0.660 

2.695 

Fig  22 

Fig  23 

0.803/CO12 

2*col4 

0.803* 
col  3 

Eqn  5 

n*Pi*a2 

1.2*n*d 

3.4  THIRD  DESIGN  EXAMPLE  (UMBRELLA  TOP-LOADED  MONOPOLE) 

The  objective  of  this  example  is  to  design  appropriate  cages  for  the  ends  of  the  top¬ 
load  radials  of  a  1200-foot  umbrella  top-loaded  monopole.  The  design  requirements  will 
be  250  kV  and  a  maximum  gradient  of  1  kV/mm.  The  diameter  of  the  wires  is  1.01 
inches. 

The  antenna  corresponds  to  VLF  Aguada  and  is  illustrated  in  figure  26.  It  has  12 
active  top-load  radials.  These  radials  operate  at  high  voltage.  Each  top-load  radial  is 
shielded  by  other  top-load  radials  and  the  tower.  Because  all  the  shielding  elements  are  at 
nearly  the  same  voltage  and  most  of  the  shielding  elements  are  the  same  size,  the  effect 
of  the  shielding  can  be  approximated  by  using  the  apparent  height  concept  introduced  in 
section  1.4.2,  Once  the  apparent  height  has  been  determined,  the  design  procedure  is  the 
same  as  for  a  single  isolated  wire. 

An  electrostatic  computer  code  (Hansen  &  Simpson,  to  be  published)  has  been  used  to 
obtain  the  charge  distribution  along  one  radial  for  a  voltage  of  100  kV  (figure  27).  The 
charge  density  is  greatest  near  the  ends  of  the  wires  where  they  are  less  shielded.  Near 
the  tower  top,  the  charge  density  is  small  because  the  radials  are  well  shielded  by  each 
other  and  by  the  tower.  Without  a  corona  ring  at  the  end  of  the  wire,  the  charge  density 
increases,  rapidly  approaching  twice  the  value  it  would  be  if  the  wire  continued  on.  VLF 
Aguada  has  large  corona  rings  at  the  wire  ends  that  mitigate  the  end  effect.  Note  that,  for 
this  case,  a  linear  charge  distribution  is  a  good  approximation. 
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Figure  25.  First  design  example,  £„/Ei  for  various  values  of  n 
(2h/a  =  28,515). 


12  TOP-LOAD  RAOIALS 


SIDE  VIEW  PUN  VIEW 

Figure  26.  Umbrella  top-loaded  monopole. 
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Figure  27.  Aguada  top-load  surfac‘d  gradient,  100  kV  with  and  without 
corona  ring. 

The  maximum  surface  electric  field  near  the  end  of  the  radials  (with  corona  rings)  for 
100  kV  is  0.696  kV/mm. 

To  determine  the  apparent  height  for  the  high-gradient  region  near  the  end  of  the 
wires,  find  the  height  for  a  single  isolated  wire  (same  diameter  and  voltage)  that  gives  the 
same  gradient.  For  a  wire  diameter  of  1.01  inches,  100  kV  and  a  gradient  of  0.696  kV/ 
mm,  the  apparent  height  is  1580  feet. 

Note  that  the  apparent  height  was  also  calculated  by  using  the  electrostatic  code  for 
wire  diameters  of  2  and  3  inches  and  found  to  be  respectively  1813  feet  and  1845  feet. 
The  apparent  height  increases  slightly  with  wire  diameter  for  this  geometrj'.  If  the  value 
1580  is  used,  the  cage  design  will  be  slightly  conservative  in  that  the  actual  gradient  will 
be  slightly  less  than  the  calculated  gradient. 

3.4.1  Given  Yfire  Size 

Design  cages  that  use  the  1.01 -inch-diameter  wire  meet  the  design  criteria. 

Using  the  value  of  given  above,  a  set  of  curves  for  versus  bla  were  calculated 
using  equation  (5)  (figure  28).  The  design  voltage  of  250  kV  has  been  used.  Note  that  in 
this  figure,  the  gradient  has  not  been  normalized. 

The  curve  for  n  =  1  is  a  special  case  for  a  single  isolated  wire  and  the  ordinate  for  this 
case  is  wire  diameter.  The  value  of  the  diameter  corresponding  to  is  determined  from 
the  value  of  the  ordinate  when  =  I  kV/mm,  and  has  the  value  1.854  inches. 
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b/4  (WIRE  DIA  WHEN  n«1) 

Figure  28.  Umbrella  top-loaded  monopole  250  kV,  1.01-inch  wires. 


There  is  no  solution  for  /i  =  2.  Note  that,  for  this  case,  the  minimum  gradient  does 
come  close  to  1  kV/mm.  Use  of  the  design  curve  of  figure  19,  with  a  =  1.01  inch/2  and 
=  1580  feet,  gives  b^/a  =  13.  The  optimum  wire  to  wire  spacing  for  the  1.01-inch  wires 
is  13  inches.  Substitution  into  equation  (5)  indicates  that  the  minimum  gradient  for  250 
kV  is  1.09  kV/mm.  Thus,  reducing  the  voltage  to  228  kV  will  satisfy  the  1  kV/mm  criteria. 

There  are  solutions  for  n  >  2,  These  solutions  are  summarized  in  table  6.  From  the 
table,  the  following  conclusions  can  be  drawn: 

(1)  The  three-wire  cage  will  satisfy  the  criteria  if  the  spacing  between  wire  centers  is 
at  least  2.1  inches. 

(2)  The  four-wire  cage  is  acceptable  if  the  center  to  center  spacing  is  at  least  1.31 
inches. 


Table  6.  Third  design  example— VLF  Aguada  with  Vtop  =  250  kV  cages  by  using 
1.0-inch-diameter  wire,  with  1  kV/mm,  apparent  height  of  wires  =  1580  feet. 


n 

a  inch 

bfa 

b  inch 

Wire  spacing 
between  centers 

S  inch 

mix 

(kV/mm) 

1 

0.927 

1.00 

0.93 

2 

0.505 

3 

0.505 

2.4 

1.212 

2.10 

1.00 

1.31 

4 

0.505 

1.84 

0.9292 

1.31 

1.00 

1.13 
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3.4.2  Minimum  Wire  Cross  Section 


The  second  design  procedure  is  used  to  find  the  cage  with  minimum  wire  cross  sec¬ 
tion,  satisfying  the  design  constraints.  The  design  curves  in  figures  22  and  23  are  used. 
The  results  are  given  in  table  7  for  250  kV. 

The  following  conclusions  follow  from  table  7. 

(1)  The  two-wire  cage  requires  wires  with  diameter  1.11  inches,  and  wire-to-wire 
spacing  of  18.7  inches. 

(2)  The  required  wire  diameter  for  n  =  3  is  0.792  inch  and  decreases  with  increasing 
n. 

(3)  For  n  s=  6,  the  required  wire  diameter  is  less  than  1/4  inch. 

(4)  The  required  cage  radius  (b^)  is  about  9.5  inches  for  all  cases. 

Table  7.  Third  design  example— minimum  wire  cross-section  wire  diameter  and 
cage  radius  for  various  values  of  n,  250  kV,  =  1580  feet,a,^£  *  0.927  inch, 
2hJa,^E  “  40,906,  £o  »  1  kV/mm. 

”  ^miJOeqE  ^  d  inch  b  inch 

1.00 
1.00 
1.00 
1.00 
1.00 
1.00 


1 

1 

0.927 

1.854 

0.927 

2 

1.67 

10.05 

0.555 

1.110 

9.316 

3.216 

3 

2.34 

10.16 

0.396 

0.792 

9.416 

4.723 

4 

3.04 

10.28 

0.305 

0.610 

9.526 

5.698 

5 

3.71 

10.32 

0.250 

0.500 

9.567 

6.367 

6 

4.39 

10.38 

0.211 

0.422 

9.622 

6.863 
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APPENDIX  A 


SINE  PRODUCT  IDENTITY 

The  sine  product  identity  to  be  proven  is  given  below, 

N-1 

n  [2sin(jbr/N)]  =  N  .  (A-l) 

k  =  l 

Euler’s  identity  for  Sin(X)  is  given  by 
Sin(X)  =  {e"'- 0/(20  . 
where  i  *  (-  1)*^*  . 

Substituting  Euler’s  identity  into  the  tefthand  side  of  equation  (A-l),  and  letting  X  = 
ht/N  gives 

N-1 

n  12{e“- 0/(20]  . 

k^l 

Operating  on  this  result  gives 
N-1 

=  (l/0^"‘  n  {e‘*-0  . 

*  =  1 


N-1 

=  (l/0^’‘  n  , 

k^l 

N-1  N-1 

=  {l/iY~\l  (-iX)}  X  n  {e^‘*-l}  , 
il:»l  k=l 


N-1  N-1 

=  (IAO  exp{S  i-ikji/N)}  X  n  <e"'*-l}  , 
Jfe*l  k=l 


N-1 

=  (-  0^'*  exp{(-  m/N)(N(N  -  l)/2)}  x  H  -  1} 

ifc=l 


N-1 

=  (-o'""’  exp{(-/;r)((N-l)/2)}x  n  {e'"-l}  . 

ik  =  l 


A-l 


N~1 

=  (exp{-/;r/2})^-'  n 

k^l 


N-1 

=  (_,-)^-»  (_,/-!  n  . 

/fe  =  l 


A/-1 

=  (-1)^"^  n  . 

ife  =  i 

*  (_  1)^-*  (- 1)^-*  ^  {1  -6^“}  . 
*  =  1 

and  finally, 

N-1 

=  n  . 

k^l 


or 

N-1 

=  n  .  (A-2) 

iS:  =  1 

The  set  of  exponentials, 

^2kx/N  forjt=lto  N-1. 

are  the  Nth  roots  of  1  (except  for  1  itself).  These  roots  are  sometimes  called  the 
cyclotomic  numbers.  The  other  Nth  root  of  1  is  1  itself.  This  can  be  seen  by  the  following: 

1  a  for  k  =  integer  . 

Taking  the  Mh  root, 

jl/N  ^  gOkZit)/N 

These  roots  are  represented  on  the  complex  plane  by  points  on  the  unit  circle  at  angles 
of  klir/N.  There  are  a  total  of  N  unique  roots  (angles),  because  as  k  increases,  the  angle 
rotates  around  to  a  unique  angle.  When  ik  =  N,  the  angle  is  2ir,  and  the  value  of  the  root  is 
1.  For  jk  >  N,  the  roots  (angles)  repeat. 

These  roots  are  solutions  of  the  equation 

-  1  «  0  . 

Using  these  roots  to  construct  the  equivalent  polynomial  gives 
2^  -  1  =  •  •  •  (Z- 1)  . 


Rearranging  terms  gives 

=  (Z^-i)/(Z-l) 


N-1 

n  (Z-e‘'*^n  =  (Z^-1)/(Z-1)  . 
jfc=  1 

Note  that  the  polynomial  Z^  -  1  can  be  factored  as  follows: 
Z^  -  1  =  (Z-1)(Z^'^  +  Z^‘^  +  •  •  •  +Z  +  1)  . 
Substituting  into  the  above  equation  gives 

^n\z-e‘'*^''^)  =  (Z^‘^  +  Z^‘^  +  •  *  •  +  Z+1)  . 
jt=l 

Setting  Z  s  1  in  the  above  equation  gives 
N-1 

n  =  N  . 

k  =  l 

Substitution  of  the  above  into  equation  (A-2)  completes  the  proof. 


A-3 


APPENDIX  B 


INFINITE  WIRE  GRID  CAPACITANCE 

The  geometry  is  an  infinite  set  of  wires  parallel  to  ground  and  each  other.  They  all 
have  the  same  radius  a  and  height  above  ground  h.  They  are  all  separated  by  the  same 
distance  S.  The  geometry  is  illustrated  in  figure  B-1.  In  the  figure,  the  slant  range  between 
the  image  of  one  wire  and  other  wires  is  denoted  by  ,  where  the  subscript  refers  to  the 
number  of  S  spacings  between  the  wires.  Z)*  is  given  by  the  following  equation. 

D,  =  7(2/1)"  +  (*5)"  . 

•  •  •!©  oooooooo*®* 


1 

1 
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1  j  /  / 

:  ;  / 

'  '  ;  /' 

WIRE  RADIUS  >  a 

1  «  ,  / 

1  .  ,  / 

'  '  '  / 

1"' 

SEPARATION  -  s 

LENGTH  <=  h 

II// 

«/// 
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Figure  B-1.  Geometry  for  an  infinite  above-ground  wire  grid. 


The  wires  all  have  the  same  voltage  V  and  by  symmeiry  they  all  have  the  same  charge 
q.  The  potential  at  one  wire  is  the  sum  of  the  potential  from  the  charge  on  the  wire  and 
the  potential  due  to  the  charge  on  all  the  other  wires.  This  potential  is  given  by  the 
following  equations: 


2;reoV 

9 


00 

+  2 


or 
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The  summation  term  can  be  rewritten  as  a  f  roduct 


The  infinite  product  in  the  above  equation  has  a  closed  form.  The  derivation  starts 
with  the  following  series  (Jolley,  1961): 


Solving  for  theta  to  make  the  products  equal, 


-dJ  ■  (fJ 

9 .  M 


Substituting  for  theta  in  equation  B-1,  gives  the  closed  form, 

.  (  .  (Inh 

j  =  1  V  \.U  J  27th  f27th\ 


(“) 


(B-1) 


Using  the  closed  form  gives  the  following  equation  for  the  capacitance  per  unit  length 
of  one  of  the  wires. 


B-2 


From  the  above  equation,  it  is  seen  that  each  wire  has  the  same  capacitance  it  would 
have  if  it  were  elevated  to  a  height  of  ■ 
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APPENDIX  C 


IMAGE-THEORY  FORMULATION 


1.  GENERAL 

The  method  of  images  will  be  used  to  calculate  the  capacitance  and  electric  fields  in 
the  vicinity  of  a  set  of  parallel  cylinders  above  a  ground  plane.  The  normal  approach  is  to 
use  a  series  of  images  to  produce  zero-potential  on  the  surface  of  two  cylinders.  This 
method  results  in  a  nonconvergent  infinite  series.  A  convergent  Image  Theory  will  be 
developed  by  using  dipoles.  This  approach  accounts  exactly  for  the  field  distortion  due  to 
the  proximity  of  the  cylinders,  but  does  not  account  for  the  field  distortion  due  to  the 
ground  plane.  It  is  essentially  exact  when  the  cylinders  are  at  a  large  distance  above 
ground  with  respect  to  the  radius  of  the  wires. 

2.  THEORY 

The  developnient  of  the  Image  Theory  starts  with  the  simple  problem  of  a  line  charge 
q  parallel  to  a  cylinder  in  free  space. 

The  geometry  for  this  problem  is  shown  in  figure  C-1.  It  can  be  shown  that  the  bound¬ 
ary  condition  (tangential  £  =  0),  will  be  satisfied  by  using  an  image  located  at  a  point 
inside  the  cylinder  boundary  having  the  same  magnitude  but  the  opposite  sign  as  the  line 
charge.  The  image  location  is  offset  from  the  center  of  the  cylinder  toward  the  line 
charge,  a  distance  given  by 

o  =  a2/D, 

where  o  is  the  offset,  a  is  the  radius  of  the  cylinder,  and  D  is  the  distance  between  the  line 
charge  and  the  center  of  the  cylinder. 

Since  the  cylinder  must  remain  at  neutral  charge,  another  charge  is  placed  at  the 
center  of  the  cylinder  having  the  same  sign  as  the  line  charge.  The  image  charge  plus  the 
charge  at  the  center  form  a  dipole  having  no  net  charge. 

The  electric  field  at  the  point  on  the  cylinder  opposite  the  line  charge  is  given  by 
E  *  {9/(2Tr£o)}  {!/(£)  +  a)  +  l/c  -  l/(a  +  n)}  . 

Substituting  for  o  gives 

£  »  {9/(2ir£o)}  {1/(D  +  a)  +  Ma  -  \l{a  +  o)} 

*  {9/(2ir£o)}  {!/(£>  +  a)  +  o/(a2  +  ao)} 

=  {qlilTTEo)}  {1/(D  +  a)  +  o/(a2  +  ao)} 

=  {q/ilTiEo)}  {2/(Z)  +  fl)}, 

which  is  exactly  twice  the  value  of  electric  field  that  would  be  at  that  position  if  the 
cylinder  were  removed. 
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Figure  C-1.  Geometry  for  line  charge  and  cylinder. 


The  series  of  images  for  the  two  cylinders,  one  charged  and  one  neutral,  is  constructed 
by  starting  with  a  charge  at  the  center  of  one  cylinder,  leaving  the  other  cylinder  neutral, 
and  constructing  an  infinite  series  of  dipoles.  The  first  dipole  is  in  the  neutral  cylinder 
and  is  the  same  as  described  above  for  the  line  charge  cylinder  combination  with  distance 
D  being  the  same  as  2b,  the  separation  between  cylinder  centers.  Each  charge  will  be 
described  by  the  sign  charge  and  offset  distance.  For  example,  the  charged  cylinder  starts 
out  with  a  +  charge  located  at  the  center  (0  offset)  designated  by  “+,  0.”  This  gives  rise  to 
two  charges  on  the  neutral  cylinder  designated  by  a2l{2b)”  and  +,  0”  (figure  C-2). 

Each  of  these  charges  gives  rise  to  two  charges  on  the  charged  cylinder.  The  charge 
+,0  gives  -,  a2/(2A),  and  +,0.  While  the  charge  -,a2l(2b)  gives  +,a2/{2b-a2l2b)  and  -,0. 
Note  that  the  two  charges  at  the  center  of  the  cylinders  (zero  offset  +.0  and  -,0)  cancel 
out.  Also  note  that  the  offset  for  the  next  image  of  an  image  is  a  simple  function  of  the 
offset  of  the  image  itself. 

From  the  description  and  the  figure,  two  important  facts  can  be  deduced:  (1)  all  image 
charges  at  the  center  of  the  charged  cylinder  will  cancel  out  and  all  image  charges  at  the 
center  of  the  neutral  cylinder  will  cancel  out  except  the  first  one;  and  (2)  the  offsets  can 
be  calculated  by  a  recursion  formula. 

If  the  first  offset  is  given  by 

01  «  a2/(2b), 

then  the  other  offsets  are  calculated  by 

Oi  =  a2/(26-o,._i). 
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Figure  C-2.  Cylinder-to-cylinder  images. 

The  distance  between  the  offsets  determines  the  dipole  moment  of  that  pair  of 
charges.  The  series  will  converge  if  each  succeeding  image  dipole  moment  is  less  than  the 
source  dipole  moment.  It  can  be  shown  that  the  offset  converges  to 

Limit  Oj  =  b  ±  ib2  -  a2)i/2. 
i  «>  infinity 

Since  the  offset  converges  to  a  fixed  value,  the  distance  between  successive  offsets 
approaches  zero,  and  the  series  converges. 

Table  Cl  illustrates  the  construction  of  the  series,  for  both  cases  with  one  cylinder 
charged,  the  other  cylinder  neutral. 


Table  Cl.  Series  of  dipoles  for  the  two-cylinder  case,  one  cylinder  charged, 
the  other  cylinder  neutral. 


Cylinder  1  (Charged) 

Cylinder  2  (Neutral) 

First  charge 

+.0 

*> 

+.0 

-.01 

=> 

-.01 

+,02 

=> 

+,02 

-.03 

=> 

-.03 

+,04 

=> 

+.04 

-.05 

=> 

Cylinder  1  (Neutral) 

+.0 

-.01 

=> 

Cylinder  2  (Charged) 
First  charge  +,0 

=> 

+,02 

-.03 

=> 

-.01 

+.02 

=> 

+.04 

-.05 

=> 

-.03 

+,04 

=> 
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From  table  Cl,  it  is  seen  that  the  sign  for  charges  at  the  even  offset  locations  is  +,  while 
the  sign  for  the  charges  at  the  odd  offset  locations  is 

When  both  cylinders  are  charged,  the  solution  is  by  superposition,  the  sum  of  the 
charges  on  each  cylinder  for  both  cases  given  in  table  Cl.  To  insure  convergence,  the 
charges  must  ’  c  taken  in  dipole  pairs.  Table  C2  illustrates  the  sum  of  the  two  series. 


Table  C2.  Sum  of  dipole  pairs:  both  cylinders  charged. 


Cylinder  1  (Charged) 

Cylinder  2  (Charged) 

+2,0 

+2,0 

-2,01 

+2,02 

-2,01 

+2,02 

-2,03 

+2,04 

-2,03 

t2,04 

-Pi 

Several  important  points  can  be  seen  by  examination  of  table  C2. 

(1)  The  series  of  charges  are  the  same  on  both  cylinders. 

(2)  The  charge  has  a  magnitude  of  2  at  every  offset  except  for  the  last  one  used.  The 
last  charge  location  has  a  magnitude  of  1,  is  always  an  odd  offset  and,  therefore,  is 
negative.  This  is  a  result  of  taking  the  charges  in  dipole  pairs. 

(3)  The  total  charge  on  each  cylinder  is  q,  the  originally  assumed  charge. 

3.  CAPACITANCE 

The  capacitance  of  the  pair  of  cylinders  is  given  by 
C  =  2q/P, 

where  q  is  the  charge  on  one  cylinder,  and  P  is  the  potential  of  the  cylinders  (assumed  to 
be  the  same  for  both  cylinders). 

The  potential  of  a  cylinder  can  be  calculated  by  summing  the  potential  at  a  point  on 
the  surface  due  to  each  charge.  The  point  used  to  calculate  the  potential  will  be  on  the 
surface  where  the  cylinders  are  closest  together.  The  potential  at  this  point,  due  to  the 
charges  on  both  cylinders  at  the  ith  offset  location,  is  given  by, 

V,  .  {2qi-l)‘/(2ltE^}M2h/(a-0,)  *  In {((2*)*  +  (2b)Y^n2l,-a-o,)]i  . 

The  total  potential  is  calculated  by  summing  the  above  terms  for  an  even  value  of  i, 
and  then  adding  one  more  term  /2.  To  simplify  interpretation  of  the  results,  the 
capacitance  is  often  normalized  to  the  capacitance  of  a  single  wire. 

4.  ELECTRIC  FIELD 


The  maximum  electric  field  will  occur  on  the  surface  of  the  cylinders  at  the  point 
where  the  cylinders  are  farthest  apart.  It  is  of  interest  to  calculate  also  the  minimum  field 


that  occurs  on  the  surface,  which  will  be  at  the  location  where  the  cylinders  are  closest 
together.  The  value  of  the  field  at  these  points,  due  to  the  charges  on  both  cylinders  at  a 
particular  offset,  is  given  by 

Maximum  Surface  Field 

£,  •  {2g(-I)7(2ir£o))  ((l/(<2  +  a,))  +  (l/(2i>  *  a  -  «,,)). 

Minimum  Surface  Field 

=  {2?(-1)'/(2it£o)}  l(l/(a  -  o.))  -  (1/(2*  -  a  -  o,)]. 

Again,  the  total  field  is  given  by  summing  the  terms  for  an  even  value  of  i,  and  adding 
one  more  term,  to  account  for  the  last  term  of  the  series. 

To  simplify  interpretation  of  the  results,  the  electric  field  is  normalized  to  the  value 
for  a  single  wire  with  radius  a.  A  simple  BASIC  computer  program  for  calculating  the 
normalized  electric  fields  and  the  normalized  capacitance  is  included  as  appendix  F. 


C-5 


APPENDIX  D 


DERIVATION  OF  COMPLETE  GEOMETRIC  THEORY 

The  derivation  for  the  complete  geometric  theory  is  the  same  as  for  the  simple  geo¬ 
metric  theory  except  that  the  fields  from  the  other  wires  are  calculated  on  the  surface  of 
the  wire  instead  of  at  the  center  of  the  wire.  This  leads  to  a  more  complicated  geometric 
factor. 

The  assumptions  used  in  this  derivation  are  summarized  below: 

1.  The  charge  on  each  wire  is  the  same  and  located  at  the  center  of  each  wire; 

2.  The  fields  from  the  images  of  the  wires  in  the  ground  plane  are  negligible;  and 

3.  The  fields  in  the  vicinity  of  the  wire  where  the  field  is  being  calculated  are 
approximately  uniform  and,  therefore,  the  effect  of  the  wire  can  be  accounted  for 
by  doubling  the  fields  calculated  as  if  the  wire  were  not  there. 

The  geometry  for  the  general  wire  cage  is  redrawn  in  figure  D-1.  For  convenience 
only,  the  first  wire  and  jth  wire  are  shown.  The  various  angles  involved  are  defined  in  the 
figure. 

The  field  normal  to  wire  1  at  the  point  opposite  the  cage  center  is  given  by 


where  pj  is  the  distance  from  the  center  of  the  jth  wire  to  the  field  calculation  point, 
and  aj  is  the  angle  between  pj  and  the  radial  from  the  cage  center  at  the  field  calculation 
point. 

The  1/a  term  represents  the  field  on  wire  1  from  the  charge  on  wire  1,  and  the  sum¬ 
mation  term  represents  the  field  on  wire  1  from  each  of  the  other  wires.  The  factor  of  2  in 
front  of  the  summation  accounts  for  the  effect  of  wire  1  on  the  fields  from  all  the  other 
wires. 

For  this  example,  the  fields  are  calculated  on  wire  1  but,  by  symmetry,  all  wires  have 
the  same  field. 
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aj  =  x/2  -  ♦y/2 

»y  =  »  -  X/2  +  ♦y/2 

Xj  =  »  ~  aj  =  */2  -  «y/2 

Figure  D-1.  General  cage  geometry  and  angle  definitions. 

1.  SIMPLE  GEOMETRIC  THEORY 

For  this  theory,  the  fields  are  calculated  at  the  center  of  the  wire,  i.e., 
pi“‘i 

(figure  D-1), 

cos(ay)  _  cos(xj) 

(pj)  dj 

Solving  for  each  factor  gives 


D-2 


or 


dj  = 

and 

cos(x,)  =  ■ 

Thus,  each  term  of  the  series  is  the  same 

cosQc,)  _  _1_ 

(dj)  ■  2b  ' 

and  the  series  sums  to  the  following  simple  closed  form 

2y±  ,  ^  . 

The  final  form  of  the  equation  is 

*  — ] 

23C€o  tfl  b  ] 

Substitution  of  the  total  cage  capacitance  C„  gives  the  result 


VCn  fl  n-r 


E  =  i] 

liteonXa  d  J 


or 


-  Hi  ^ 

This  equation  is  called  the  simple  geometric  theory  (SGT). 


2.  COMPLETE  GEOMETRIC  THEORY 


For  the  complete  geometric  theory,  the  fields  are  calculated  at  the  surface  of  wire  1 , 


i.e.. 


Pj  =  Sj 
Otj  =  tjf;. 

The  summation  becomes 
2  ^  cosTO 
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By  using  the  definitions  of  figure  D-1  and  the  law  of  Cosines,  we  can  solve  for  Sj 

Sj^  ^  dj^  +  -  2adjCos(^j)  . 

Substitution  of  the  following 
cos(ipj)  =  ^0^(7 

dj  =  2f?sin^^j 


gives 


Sj^  =  4d^cos|-;^J  +  +  4adsin' 


2  • 


which  reduces  to 


Sj  ^  +  4b(b  +  a)  • 

Similarly,  solving  for  cos(i|rj) 
dj^  =  Sj^  •¥  -  laSjCOsitpp 


cositpj)  = 


Sj  ^  -  dj 


The  summation  terms  are 


cos(V';)  2a' 


,  -  .  ...  a  +  26sin*r^l 

+  2adj  cosjipj)  _ y2  J 


cos(ipj) 


a  +  26  sin^l 


a^  +  (46)(6  +  a)  sin 


© 


For  n  =  2  theta  sub  j  over  2  =  pi,  and  the  equation  for  maximum  surface  field  becomes 


1  If— -1— 
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For  «  =  3  theta  sub  j  over  2  =  pi/3  and  2  pi/,  and  the  equation  for  maximum  surface 
field  is 

E  =  -fl  +  1  ][  1  1 

"  “  3\fl  fl2  +  3fl6+362jLln(2Vfl')J 

For  n  -  4  theta  sub  j  over  2  =  pi/4,  pi/2,  and  3  pi/4,  and  the  equation  for  maximum 
surface  field  is 

p  _  1  ir  1  1 

~  4  la  ^ +  lab  2b'^]Y\n{2h/a')\ 

2 

3.  MODIFIED  SIMPLE  GEOMETRIC  THEORY 

The  terms  of  the  summation  can  be  rewritten  in  the  following  form: 


These  terms  are  all  of  the  form 


1 

6  +  I  +  CF(9) 

where  CF  is  a  correction  term.  If  we  ignore  this  correction  term,  the  series  has  a  simple 
closed  form,  and  gives  the  following  equation  for  the  maximum  surface  field 

t  [b(^] 

2 

This  equation  is  very  similar  to  the  simple  geometric  theory  result  and  will  be  called  the 
modified  simple  geometric  theory  (MSGT). 
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DERIVATION  OF  CHARGE  OFFSET  FORMULAS 

The  basis  of  this  derivation  is  to  assume  the  charge  on  each  wire  can  be  represented 
by  a  filament  of  charge  located  near,  but  offset  from  the  center  of  the  wire.  The  offset  is 
in  the  direction  of  a  line  between  the  centers  of  the  wires.  The  geometry  is  illustrated  in 
figure  E-1.  The  offset  is  determined  by  the  requirement  that  the  potentials  calculated  at 
two  points  on  the  surface  of  the  wires  are  equal.  The  points  are  located  where  the  exten¬ 
sion  of  the  line  between  the  centers  of  the  wires  intersect  the  surface  of  the  wire.  These 
are  the  points  where  the  wires  are  closest  together  and  farthest  apart.  The  offset  is 
assumed  to  be  in  a  direction  such  that  the  charge  locations  move  towards  the  outside  of 
the  wires. 
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Figure  E-1.  Geometry  for  charge  offset  derivation. 


The  equations  for  the  potential  at  these  points  are 
Pi  «  {9/(2itEo)}  [ln(2A/(a  +  o))  +  ln(2fi/(2b  +  o  -  a))] 

Pz  =  {9/(2ir£o)}  [ln(2A/(a  -  o))  +  ln(2/i/(2f>  -k?  +  o))]. 

Setting  these  potentials  equal  leads  to  the  following  equation: 

\n{2hl[a  -f  o))  +  ln(2A/(2f)  +  o  -a))  =  \n{2h/{a  -  o))  +  \n(2h/(2b  +  o  +  a)). 

Eliminating  the  logarithms  and  reducing  as  follows  results  in  a  quadratic  equation 
for  0. 
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{2h)2l[{a  +  o)  (26  +  0  -  a)]  =  {2h)2l[{a  -  o)  (26  +  o  +  a)] 

(a  +  o)  (26  +  0  -  a)  =  (a  -  o)  (26  +  a  +  a) 

2a6  +  ao  -  a2  +  2bo  +  o2  -  ao  =  2ab  +  ao  +  a2  -  2bo  -  o2  -  ao 
2o2  +  46o  -2a2  =  0 

a2  +  26a  -a2  =  0.  (E-1) 

The  quadratic  equation  for  o  has  the  solutions 
a  =  -6  ±  (62  +  a2)i/2. 

The  capacitance  of  the  pair  of  wires  is  given  by 
C  =2^/P. 

whete  P  is  the  potential,  either  P\  or  P2  since  by  definition  they  are  equal. 

Substituting  the  value  for  the  offset  a  into  the  equation  for  either  or  P2  gives  the 
following, 

P  =  {q/{2TtEo)}  [ln((26)2/(a2  +  26a  -  a2  +  26a)]  . 

Substituting  equation  (£-1)  gives 
P  =  {q/(2TTEQ)}  Iln((26)2/(26a)]  , 

which  gives  the  same  potential  as  if  the  calculations  were  done  with  the  charges  located  at 
the  center  of  the  wires.  The  capacitance  for  the  two-wire  case  is  the  same  as  the  SGT 
given  by 

C  =  47rEo/[ln((2h)2/(2ba)]  , 

»  2ir£o/(ln((2/i)/(26a)i/2]  . 

The  maximum  surface  electric  field,  located  at  point  2  (figure  E-1)  is  calculated  by 
using  Gauss’  law  and  the  value  of  the  offset  as  follows, 

E  =  {9/(2ir£o)}  {l/(a  -  a)  +  1/(26  +  a  +  a)}  . 

Reducing  the  above  equation  gives 

£  =  {9/(27r£o)}  {(26  +  a  +  0  +  0  -  o)/((o  -  a)  (26  +  a  +  o))}  , 

=  {g/(2Tr£o)}  {2(6  +  a)/(2a6  +  a2  +  ao  -  26o  -  ao  -  o2)}  , 

=  {9/(2ir£o)}  {2(6  +  a)/(2a6  +  a2  -  26o  -  o2)}. 

Substitution  of  equation  (E-1)  further  reduces  this  to 
£  =  {^/(2Tr£o)}  {2(6  +  a)/(2a6)}, 

=  {g/(2ir£o)}  {(6  +  a)/(a6)}  , 

«  {g/(2Tr£o)}  {(1/a  +  1/6)}  . 

£  =  {(l/(a  -  0)  +  1/(26  +  a  +  o)/n}  [l/(ln(2/i/a0)}  • 

Finally,  substituting  for  q  gives 
£  =  V{(l/a  +  l/6)/2}  [l/ln(2/i/a')]  , 
which  is  exactly  the  same  as  the  SGT  for  the  two-wire  case. 

Note  that  dependence  on  the  offset  “0”  explicitly  cancels  out  in  the  formulas  for  both 
surface  electric  field  and  capacitance. 
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APPENDIX  F 


COMPUTER  PROGRAM  TO  FIND  AND  a  j, 

IN  BASIC 


F-1 


DECLARE 


100 


200 


300 


SUB  equlv  (h21,  *01,  aOl) 

CLS 

OPEN  "cage2.dat "  FOR  OUTPUT  AS  #i 
n  *  2 
b  •  1 

PRINT  "Location  of  Kin:muia  Gradiant  and  Nornallzed  Value" 

PRINT  "n*";  n 

PRINT  ■2h";  "  b";  "  Em";  "  aeqe";  "  aaqc";  " 

i  -  100 
h2  ■>  0 

DO  WHILE  h2  <  100000 
h2  <■  b2  -f  i 

•0  «  100 

aeqc  «  (n  *  b  ■*  (n  -  1))  “  (1  /  n) 

eon  -  (1  +  (n  -  1)  /  (b  ♦  .5))  "  L0G{h2)  /  n  /  LOC(h2  /  aeqc) 

IF  emn  <  eo  THEM 
eO  >  ean 
bO  >  b 
aeqcO  ■>  aeqc 
b  •=  b  +  1 
GOTO  100 


ELSE 

b  -  bO  -  1 
eO  «■  100 
END  IF 

aeqc  -  (n  *  b  (n  -  1) )  *  (1  /  n) 

emn  -  (1  +  (n  -  1)  /  (b  .5) )  *  LOC(h2)  /  n  /  LOG(h2  /  aeqc) 
IF  enn  <  eO  THEN 
eO  -  ean 
bO  -  b 
aeqcO  >  aeqc 
b  -  b  +  .1 
GOTO  200 


ELSE 

b  -  bO  -  .1 
eO  -  100 
END  IF 

aeqc  -  (n  ♦  b  “  (n  -  1))  -  (1  /  n) 

emn  -  (1  +  (n  -  1)  /  (b  +  .5) )  *  LOG(h2)  /  n  /  LOG(h2  /  aeqc) 
IF  eO  >>  emn  THEN 
eO  >  ean 
bO  -  b 
aeqcO  ~  aeqc 
b  -  b  +  .01 
GOTO  300 

ELSE 
END  IF 

equiv  h2,  eO,  aO 

PRINT  h2;  bO;  eO;  aO;  aeqc;  bO  /  aO;  bO  /  aeqc 
PRINT  #1,  h2;  bO;  eO;  aO;  aeqc;  bO  /  aO;  bO  /  aeqc 

IF  h2  -  1000  THEN  1  -  1000  ELSE 
IF  h2  -  10000  THEN  1  «  10000  ELSE 
LOOP 


F-2 


END 


STOP 


SUB  aquiv  (h2,  eO,  aO) 

a  ■  1 
ea  ~  eO 

'  PRINT  a,  h2,  eO 

1000  eal  -  I.OG(h2)  /  a  /  LOG(h2  /  a) 
IF  eal  >  eO  THEN 
a  «  a  +  .001 
GOTO  1000 

ELSE 

aO  >  a 

END  IF 

'  PRINT  eal,  a,  eO 

END  SUB 
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